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PREFACE 


The work described in this report was performed by the Propulsion 
Division of the Jet Propulsion Laboratory. 
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ABSTRACT 


Recent test* in the JPL Molsink facility have provided the first quanti- 
tative evidence of gaj; flows in the far upstream region of small nozzles with 
large boundary layer flow. Gas mass fluxes were measured using quartz 
crystal microbalances. Both nitrogen and carbon dioxide gases were used 
as test gases. Gas deposition rates on the order of 100 monolayers per 
minute were detected 13 inches upstream of the nozzle exit plane. It is sig- 
nificant to note that the crystals detected gases considerably beyond the 
inviscid Prandtl-Meyer turning angle. Tests are currently continuing using 
improved cryogenic quartz crystal systems and additional types of gases. 

The data from these tests will be essential in the formulation of scaling laws 
and analytical prediction methods for viscous plume behavior. 


JPL Technical Memorandum 33- 540 


vii 



INTRODUCTION 


The interaction of plumes with spacecraft components can cause prob- 
lems in such diverse areas as flight mechanics, materials and structures, 
communications, and scientific experiments. The possible problems are 
aggravated if a large percentage of the plume is composed of nozzle boundary 
layer gases. Because these have a lower velocity than the nozzle core gases, 
they have a higher potential of expansion and will cause a back, or upstream 
flow, of the exhaust gases. Although analytical methods, such, as the method 
cf characteristics, are available to predict plume characteristics for nozzles 
with negligible boundary layers, no equivalent proven methods are available 
for those with large boundary layers. Because the latter type of nozzles are 
typical of those used aboard unmanned spacecraft today, an experimental and 
analytical program to predict the nature of these plumes has been undertaken. 
The program has been described and some initial results have been reported 
in Reference 1. Recent tests using nitrogen and carbon dioxide gas have 
resulted in the first data on the amount of gas in the backflow region. Tests 
are continuing and the data reported in this TM, correlated with that now 
being generated, will result in the formulation of scaling laws and analytical 
techniques which will predict the distribution of gases in viscous plumes. 
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FACILITY AND TEST APPARATUS 


Facility Characteristics 

The Molsink is an ultra- high vacuum facility which consists of three 
concentric chambers. The innermost chamber is a sphere approximately 
3. 05 m (10 ft) in diameter, maintained at a temperature between 15 and 
20 °K with gaseous helium. The walls are wedge-shaped, resembling an 
anechoic chamber, with a total surface area of approximately 186 m2 
(2000 ft^). The chamber walls are also coated with titanium, which acts as 
“getter” material to trap helium and hydrogen that are not cryopumped by 
the 1 5°K surfaces. The facility is described in more detail in References 2, 
3, and 4. A cross section of the Molsink is shown in Figure 1. 

Quartz Crystal Microbalances (QCM) 

Quartz Crystal Microbalances were positioned at various locations 
along the Molsink walls and used as gas detection devices. The QCM meas- 
ures the mass per unit area adsorbed or desorbed on a piezoelectric quartz 
crystal surface through a change in the resonant frequency of the crystal. 

For a given crystal cut angle, there is a range of temperature for which the 
frequency change depends linearly on the rate of mass deposition. The 
QCM systems were orientated with their sensing surfaces perpendicular to a 
radial line drawn from the nozzle lip. Each QCM system included four 
crystal sensing surfaces, 1. 3 cm in diameter, with their centers approxi- 
mately 1. 6 cm apart. Figure 2 is a photograph of the QCM system sensing 
heads as installed in the Molsink. 

.The QCM systems are driven by oscillators. Not all of the installed 
systems operated satisfactorily because of crystal/oscillator coupling prob- 
lems. To condense and trap nearly all of the impinging gas the crystal sur- 
faces must be maintained at low temperatures. The temperatures required to 
condense various gases are indicated in Table 1, which is taken from Refer- 
ence 5. For example, to trap nitrogen at a pressure of 1. 33 X lO"*^ N/m^ 
(10“^ torr), the crystal temperature must be below 29°K. Net all of the 
crystals were cooled to temperatures low enough to condense nitrogen. 
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The experience gained in these tests has been used to design and con- 
struct tmproved QCM systems which are operating successfully at low 
(<Z9°K) temperatures during present tests. The relative location of the 
test nozzle and the QCM systems from which data was obtained is shown 
in Figure 3. 
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TEST RESULTS 


Tests were conducted with both nitrogen and carbon, dioxide. The gas 
was released into the Molsink through a conical nozzle while the frequencies 
of the quartz crystals were recorded approximately every 30 seconds. The 
frequencies of the quartz crystals were also recorded when no gas was flow- 
ing in order to determine the mass of gas, previously deposited, leaving the 
crystal surface. 

Test data was obtained for carbon dioxide gas at two locations, 0. 33 m 
(13 in.) upstream of the exit plane (163 deg from the nozzle axis) and near 
the exit plane approximately 0.76 m {30 in.) from the nozzle. Because the 
crystal sensing surface of the QCM system upstream from the nozzle was 
not cold enough to condense the nitrogen, data for was obtained only near 
the exit plane. 

The test results, in terms of change of mass on the sensing crystal 
surfaces, as a function of time are shown in Figures 4, 5, and 6 and are 
tabulated in Table Z. Figures 7 and 8 show the change of mass on the crystal 
as indicated by the QCM system during the gas off portion of the test and 
confirm the fact that the only significant change in frequency occurred during 
the time that gas was actually flowing through the nozzle. This verified the 
fact that desorption of the condensed gases from the crystal surfaces, during 
the test, was negligible. 

Test results also indicate that the amount of gas found in the backflow 
region is strongly related to the proportion of the total nozzle gases in the 
boundary layer. In isentropic flow, the mass flux at a fixed location in the 
plume, can be assumed to be directly proportional to the nozzle plenum 
pressure. The test data, however, indicates that for a ratio of plenum 
pressures of 4. 5:1, the mass flux ratio in the backscatter region is only 
2. 3:1. An integration of the experimentally determined Mach number pro- 
file at the exit plane of the nozzle shows that for a plenum pressure of 
2. 76 X 10 4 N/m^ (4 psia) the boundary layer flow is approximately 32 per- 
centage of the total mass flow, while at a pressure of 1 . 24 X 10 ^ N/m^ 

(18 psia) it is only 16 percent of the flow. Therefore, it would seem that the 
quantity of gas in the backscatter region is not directly proportional to the 
nozzle total mass flow and that the reduction in the boundary layer at the 
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higher plenum pressures and Reynolds numbers results in a smaller per- 
centage of the nozzle total gas flow expanding into the backflow region. 
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DATA ANALYSIS 


Data Reduction Procedure 

The raw data in the form of quartz crystal frequencies is converted to 
mass per unit area changes on the crystal surface by means of the following 
theoretical formula: 


Arn 


4. 35 x IQ 
? 

F~ 


(AF) 


Am = change in mass flux, gm/cm 

F = resonant frequency of crystal, MHz 
AF = change in frequency, Hz 

The key assumption in the above formula is that the temperature of the 
crystal remains unchanged. The temperature sensitivity is a function of the 
angle at which the crystal is cut, and the crystals used in the QCM systems 
from which data was obtained are relatively insensitive to temperature 
changes. Therefore, it was not necessary to make temperature corrections. 

Nozzle Characteristics 

The flow characteristics of the nozzle used to generate the plume were 
measured during tests in the JPL space simulator. The nozzle is conical, 
with a 60:1 area ratio and 15° wall half angle. The percentage of the exit 
area occupied by the boundary layer increases rapidly with decreasing 
chamber pressure, as determined experimentally and indicated in Figure 9. 

It is clearly evident that viscous effects cannot be neglected in predicting tbe 
plume gas distribution. 

Plume Maps 

Inviscid plume gas distributions can be calculated quite accurately using 
the method- of- character! sties (MOC) which is valid for supersonic flow. 

The plume maps shown in Figures 10 and 11 were calculated for nitrogen 
and carbon dioxide. Also indicated on the figures is the Prandtl-Meyer 
turning angle based upon the nozzle exit Mach number, as calculated by the 
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MOC program. This is the maximum theoretical turning angle for the flow 
(or the angle at which the calculated density is zero) from this nozzle neglect- 
ing boundary layer effects. 

The mass flux data from the quartz crystals can be translated into 
density if it is assumed that in the region of the quartz crystals the mole- 
cules have reached their maximum velocity (the velocity at which the total 
enthalpy has been completely converted to kinetic energy). The experimental 
data, thus converted, has also been plotted in Figures 10 and 11. 

In Figures 12 and 13 the mass flux data is presented directly, thereby 
eliminating any error that may be due to deviations from adiabatic conditions 
which were assumed in converting the mass flux data to equivalent densities 
in Figures 10 and 11. The inviscid plume data shown is non- dimen sionalized 
and indicates the mass flux off the centerline as a function of the centerline 
mass flux in the plume far field where the QCM systems are located. 

The juxtaposition of the experimental data and the analytically derived 
inviscid plume data, in Figures 10, 11, 12 and 13, emphasizes the signi- 
ficant errors that can occur in plume behavior predictions when viscosity 
effects are neglected. The flow indicated by the QCM system, even though 
small, represents a significant percentage of the nozzle exhaust gases. A 
conservative calculation, assuming the average mass flux detected by the 
QCM systems is constant everywhere in the backflow region, at equal dis- 
tances from the nozzle exit, shows that this flux represents 1. 5% of the 
nozzle total mass flow in case of the carbon dioxide tests. For purposes of 
the calculation, the backflow region was defined to be the region upstream of 
the Prandtl-Meyer turning angle. 
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CONCLUSIONS 


The significance of these t*- sts is that they offer the first data of the 
amount of gas in the far upstream, or backflow, region of small nozzles. 
Despite the fact that this initial data is too limited to allow the formulation 
of accurate scaling laws , a number of conclusions can be made: 

(1) The amount of gas in the region outside the Prandtl-Meye r angle 
is large enough to affect sensitive coatings on spacecraft instru- 
mentation. In terms of monolayers, for example, the amount of 
carbon dioxide gas deposited on the #1 QCM system was on the 
order of 100 monolayers pe'r minute. The small amount of gas 
required to coat a surface with several monolayers can be seen 
by referring to Table 3 which is a tabulation of the molecular 
properties of several gases from References 6 and 7. The mean 
free path of a molecule is proportional to the absolute tempera- 
ture and inversely proportional to the pressure. Table 3, which 
also lists the free mean path of the test gases, indicates that the 
Knudsen number (based on the sensing element diameter of 
appr oximately 1. 3 cm) in the vicinity of the QCM systems is 
large enough for free molecular flow conditions to exist. 

(Z) The proportion of the nozzle gases in the backflow region is a 
strong iunction of the nozzle boundary layer thickness. An in- 
crease in the nozzle gas flow of 4. 5 coupled with a decrease in 
the percent of gases in the boundary layer of 50% resulted in an 
increase of Z. 3 in the gases in the backflow region. The amount 
of gas in the backscatter region can therefore be reduced by 
nozzles, which by the nature of their design, inhibit the growth 
of boundary layers. This is one of the items being studied in an 
analytical effort paralleling the present experimental work and 
whose primary objective is to analytically predict the properties 
of viscous plumes. 

(3) Finally, it can be stated that these tests have provided initial 
data which will help the formulation ci analytical techniques 
required to accurately predict the characteristics of 
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viscous plumes. This is necessary to insure the compatibility 
of the spacecraft, its propulsion system, and its science package 
so that the mission objectives are not inadvertently compromised. 
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Table 2. Rate of mass deposition of CO^ and N^ plume gases 


Gas (1) 

O 

o 

rj 

C °2 

N 2 

N z 

Nozzle Plenum Pressure, 
N/m^ (psia) 

2. 8 x 10 4 
(4) 

2 . 8 x 10 4 
(4) 

2. c X 10 4 
(4) 

1. 2 x lo’ 
(181 

Quartz Crystal Location^ 

#1 

# 2 

=z 

- 2 

Mass Flux, kg/m“ ~ s 

1. 1 X 10" 6 


2. 7 X 10' 6 

1 . 1 x 10'° 

2. 5 X 10"° 


Note: 


( 1 ) 


Gas injected through 15° 
diameter of Z. 54 X 10"^ 


conical nozzle with 60:1 
m (0. 100 in. ) . 


area ratio and a throat 


'‘"'Quartz Crystal location shown in Figure 2. 

rr 1 Approximately 0. 33 m (13 in.) upstream of exit plane. 

=T 2 Approximately at exit plane and 0. 76 m (30 in. ) from nozzle 
centerline . 


Table 3. Molecular properties of gases 




C °2 

N 2 

H 2 

1 . 

. 24 

Mass of molecule, 10 g 

73 

46. 5 

3. 35 

z. 

Numb, r of molecules per 
cm^ in a monolayer, 10^ 

5. 34 

8. 10 

15. 22 

3. 

Mean free path at 15°C and 
0. 133 N/xn* (1 p Hg), cm 

3. 18 

4.77 

8. 97 

4. 

2 

Mass of monolayer per cm , 
10* 8 g 

3.9 

3.8 

0. 5 

i 
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Fig. 2. QCM System Sensing Heads, as Installed in the Molsink 
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Fig, 5 , Quartz Crystal Mass Deposition- 
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Fig. 6. Quartz Crystal Mass Deposition- 
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Flow Terminated 
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NOZZLE PLENUM PRESSURE, N/m 2 (psia) 


Fig. 9. Nozzle Boundary Layer Growth 
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Carbon Dioxide Plume Map 





Plume Boundary Data--CO 
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